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Abstract 

Using a sample of 58 million J /ip events recorded in the BESII detector, the decay 
J/ip — ► ujtt + tt^ is studied. There are conspicuous o>/2(1270) and &i(1235)7r signals. 
At low 7T7T mass, a large broad peak due to the a is observed, and its pole position 
is determined to be (541 ± 39) - i (252 ± 42) MeV from the mean of six analyses. 
The errors are dominated by the systematic errors. 



PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 



1 Introduction 



Two important aspects of strong interaction physics at low energies are color 
confinement and spontaneous breaking of chiral symmetry. Better experimen- 
tal knowledge on the ++ low energy hadron spectrum is very important to 
the understanding of QCD in the nonperturbative region. For example, the 
determination of the lowest lying ++ hadron, named the /o(600) or a, will be 
very helpful in understanding how QCD realizes chiral symmetry. 

There has been evidence for a low mass pole in the early DM2 [1] and BESI 
[2] data on J /if) — > ujtt + 7t~ . A huge event concentration in the 1 = S-wave 
mi channel was seen in the region of m nn around 500-600 MeV in a pp central 
production experiment [3]. This peak is too large to be explained as back- 
ground [4] . There have been many studies on the possible resonance structure 
in 7T7T elastic scattering. A summary of these studies up to 1999 is given by 
Markushin and Locher [5]. It was later proved that the a resonance is unavoid- 
able in chiral perturbation theory in order to explain the tttt scattering phase 
shift data [6]. An analysis based on chiral symmetry and the Roy equations 
has been made by Colangelo, Gasser and Leutwyler [7] ; a light and broad res- 
onance was found with the result M - iV/2 = (470 ± 30) - i(295 ± 20) MeV 
for the pole position. Renewed experimental interest arises from E791 data on 
D+ -»■ 7T + 7r-7r+ [8]; they find M = 478i| ± 17 MeV, T = 324+^ ± 21 MeV. 

Results on J /if) — ► c<J7r + 7r~ from 5.8 x 10 7 J/if) events collected with the up- 
graded BES (BESII) detector are presented in this paper. The upgrade of the 
BES detector included a new main drift chamber (MDC) and a new time- 
of-flight (TOF) system. A detailed description of the BESII detector is given 
in Ref. [9]. It has a cylindrical geometry around the beam axis. Trajectories 
of charged particles are measured in a vertex chamber (VC) and the main 
drift chamber (MDC); these are surrounded by a solenoidal magnet provid- 
ing a field of 0.4T. Photons are detected in a lead-gas Barrel Shower Counter 
(BSC). Particle identification is accomplished using time-of-flight (TOF) in- 
formation from the TOF scintillator array located immediately outside the 
MDC and the dE/dx information in the MDC. 

Partial wave analyses (PWA) are performed on this channel using two meth- 
ods. In the first method, the whole mass region of M n + n - which recoils against 
the u) is analyzed, the u) decay information is used, and the background 
is subtracted by sideband estimation. For the second method, the region 
M n + n - < 1.5 GeV is analyzed, and the background is fitted by 5n phase space. 
In both methods, different parametrizations of the a pole are also studied. 



2 Event selection 



Here, the selection of ujtt + 7t~ events is outlined, where the u is observed in 
its 7t + 7t~7t° decay mode. Candidate tracks are required to have a good track 
fit with the point of closest approach of the track to the beam axis being 
within the interaction region of 2 cm in ^Jx 2 + y 2 and ±20 cm in Z (the beam 
direction), polar angles 9 satisfying |cos#| < 0.80, and transverse momenta 
> 60 MeV/c. Photons are required to be isolated from charged tracks and 
to come from the interaction point. Any photon with deposited energy lower 
than 30 MeV in the BSC is rejected. Events are required to have four good 
charged tracks with total charge zero and more than one good photon. The 
TOF and dE/dx information are used to identify pions; they largely reject 
kaons from background reactions such as K + K~ n + n~ it . 

A four constraint (4C) kinematic fit is applied under the 7r + 7r~7r + 7r~77 hypoth- 
esis, and xlc < 40 is required. Events with a 27 invariant mass \M 71 — M n o\ < 
40 MeV/c 2 are fitted with a 5C kinematic fit to 7t + 7t~7t + 7t~tt with the two 
photons being constrained to the 7r° mass. Events with \\c < 40 are selected. 
The resulting tt + tt~it mass distribution is shown in Fig. 1(a). If there is more 
than one mass combination satisfying requirements, the one closest to the u 
mass is plotted. The u signal is selected by requiring \M v +^- v o — M u \ < 40 
MeV/c 2 . Separation of ujtttt from uKK is very clean, since the kinematics of 
these two processes differ strongly. 

Fig. 1(b) shows the 7r + 7r~ invariant mass spectrum which recoils against the 
ou, and Fig. 1(c) shows the un invariant mass. The Dalitz plot of this channel 
is shown in Fig. 1(d). There is a large ^(1270) peak in Fig. 1(b) and a strong 
6i(1235) peak in Fig. 1(c). At low ttti masses in Fig. 1(b), a broad enhancement 
which is due to the a pole is clearly seen. This peak is evident as a strong 
band along the upper right-hand edge of the Dalitz plot in Fig. 1(d). Events 
produced according to phase space will be broadly distributed over the whole 
Dalitz plot region. Therefore this band does not correspond to phase space. A 
detailed Monte Carlo simulation is performed to study the background. The 
main background comes from 5-7T events without an uo and is determined to be 
(14.4 ±1.5)% using a sideband estimation. 



3 Partial wave analysis 



PWA analyses are performed. Two methods and different parametrizations of 
the a pole are used. 

For the first analysis, Method I, the whole region of 7r + 7r~ mass recoiling 
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Fig. 1. (a.) Distribution of tt + tt tt° mass, (b.) Distribution of the 7r + 7r invariant 
mass recoiling against the to. (c.) Distribution of coir invariant mass, (d.) Dalitz plot. 

against the u> is considered. The channels fitted to the data are: 



J /ip-> uj / 2 (1270) 

-u;/o(980) 
^6i(1235)tt 
^p(1450)tt 
->/ 2 (1565)w 

^/ 2 (2240)u;. 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 



Amplitudes are fitted to relativistic tensor expressions which are documented 
in Ref. [10]. For spin in tttt, two transitions from J/ijj are allowed with orbital 



angular momenta £ = and 2 in the production process. For spin 2, there are 
five amplitudes: one with £ = 0, three with £ = 2 and one with £ = 4. In 
fitting these, Blatt-Weisskopf centrifugal barrier factors are included with a 
radius of 0.8 fm. Each amplitude is fitted with a complex coupling constant 
G, following the usual isobar model. 

The subtraction of background is made using sidebands 80 MeV/c 2 wide, 
centered at M(ii + tVtt ) = 622 and 942 MeV/c 2 . Signal events in the u mass 
band are given positive weight in log likelihood and sideband events negative 
weight; the sideband events (suitably weighted by 7r + 7r~7r° phase space) then 
effectively cancel background in the data sample. Results are stable when the 
position and width of the sidebands are varied. 

In the amplitude analysis, information from the io — ► 7r + 7r~7r° decay is included 
in the tensor expressions. The u polarization is given by u a = € a p y $(k + )t \k~)' 1 (k°) s ', 
where k + , k~ and k° are the 4-momenta of n + , n~ and 7r°, and e is the fully 
antisymmetric unit tensor. Non-relativistically, this reduces to the vector prod- 
uct fc + x fc~, where k + and k~ are the three- momenta of the 7r + and the 7r~ 
in the rest-frame of the u. The spin of the u therefore lies along the normal 
to the u — ► 37r decay plane. Angular correlations between this normal and (a) 
the beam direction, (b) the 7t + tt~ pair provide a delicate separation between 
spins and 2 for the pion pair; they also identify different L values in the 
production process. 

In the second analysis, Method II, the 7r + 7r~ pairs recoiling against the u with 
M n + n - < 1.5 GeV/c 2 are studied. The covariant helicity coupling amplitude 
method is used [11] to construct the amplitudes. The u decay information is 
not included in the analysis. In the 7r + 7r~ invariant mass spectrum, the pro- 
cesses of uia, w/2(1270), and w/o(980) are included in the fit, and &i(1235)7r 
is considered in the urn invariant mass spectrum. The background is approxi- 
mated by a non-interfering 5n phase space. 

In both analyses, the /o(980) is parametrized by the Flatte formula which is 

written as: 

1 
/ 



M 2 - s- i(gip^(s) + g 2 pKK(s)) ' 



where p(s) = 2k/ y/s and k is the center of mass momentum of the n or K in 
the resonance rest frame; we take M = 0.970 ± 0.007 GeV, g x = 0.138 ± 0.010 
GeV, and g 2 /gi = 4.45 ±0.25. The full- width at half maximum is 24 ±3 MeV. 

For the ^(1270) , only two of the three possible £ = 2 amplitudes are signifi- 
cant, and the £ = 4 amplitude is small. In Method I, the ^(1270) is optimized 
as a Breit-Wigner with a width proportional to k 5 B 2 , where B 2 is the L = 2 



centrifugal barrier, while in Method II, a constant width Breit-Wigner is used 
for the / 2 (1270). The &i(1235) mass optimizes at 1225 ±4 MeV and 1231 ±12 
MeV in Methods I and II, respectively. These are consistent with the mean 
mass quoted by the Particle Data Group (PDG) [12]. The / 2 (1270) mass is set 
to the PDG value in the final fits. The ratio of D and S-wave decay amplitudes 
of the &i(1235) is consistent with the PDG value of 0.29 and is fixed to this 
value. 

We come now to the a pole. Several parametrizations of the a pole are used 
in the analyses. The first is a Breit-Wigner with a constant width: 



""- W-.-iMT^ (9 » 



Secondly, the form introduced by Zou and Bugg [13] in fitting im elastic 
scattering data, consistent with Cern-Munich data on nix elastic scattering 
[14] and with K e ^ data [15], is used for the a amplitude. 



f M*-s-iMT m (sy (10) 

" (S) " 5l ^) +S2 ^)' (11) 

9i = m^f^- 2 eM-(s-M^)/a}. (12) 



Here p nn is the usual nn phase space 2k/y/s, and k is the momentum in the irir 
rest frame. The form includes explicitly into T(s) the Adler zero at s = m^/2; 
the exponential factor cuts off the width at large s. A revised fit to tctt elastic 
data and K e ± using this formula is presented in [16]. 

In Eqn. (12), f(s) — b\ + b 2 s, where b\ and 6 2 are adjusted to reproduce the 
scattering length and effective range for im elastic scattering; these are from 
recent K e ^ data of Pislak et al. [15]. In the second term of Eqn. (11), 47r phase 



space p4,r(s) is approximated by J(l — 16/i 2 /s)/[l + exp(2.8 — s)/3.5], with 
s in GeV 2 . In practice, the 4ir width is significant only at masses above 1200 
MeV and has no effect on the a pole. 

Thirdly, we have also tried fitting the a with: 



BW - K -.-wy TJs) ^ M ^ 1 ' (13) 



which was used by the E791 collaboration [8]. For the Breit-Wigner amplitude 
with an s-dependent width shown in Eqn. (13), the pole position M does not 
coincide with the mass M where the phase shift goes through 90° [16]. 

Fourth, the form given in Ref. [17], which removes the spurious singularity 
hidden in Eq. (13), is used: 



BW - = mi -.-iV!r.W r -W = «Vi-<- < 14 > 

The fits are made using the maximum likelihood method. 

4 Results and discussion 

4.1 Method I 



In method I, when using Eqns. (10) - (12) to describe the a pole (Fit A), 
the optimum fit is obtained with M = 0.9264, g 2 = 0.0024, a = 1.082, b t = 
0.5843, 62 = 1-6663 (all in units of GeV). Fig. 2 shows the projection of the 
fit compared with data. 

The BES data determine the mass and width of the pole well, but with a 
rather strong correlation between them. Fits have been made with 46 variants 
for f(s) in Eqn. (12). Optima lie roughly along a line from M = 500 — i270 
MeV to 600 — il95 MeV: the fitted width goes down as the mass goes up. 
The optimum is at (542 ± 7 (stat) ±15 (sys) ±30 (extrap)) - i(249± 15 (stat) 
±20 (sys) ±30 (extrap)) MeV. Systematic errors arise roughly equally from (i) 
varying the choice of sidebins, (ii) varying the magnitude of the background 
under the 10 peak of Fig. 1(a), and (iii) varying the choice of small components 
in the fit (discussed below). The last error accounts for systematic errors in 
the extrapolation to the pole. This has been estimated from the 46 variants 
used to fit BESII data and therefore covers the error due to different formulae. 
The systematic errors dominate. 

The mass and width of the / 2 (1270) optimize at 1271 ± 5 MeV and 174 ± 10 
MeV, respectively. The fitted width of the b x (1235) is 195 ± 20 MeV; this is 
distinctly larger than the PDG value of 142 ± 9 MeV. A fit with the PDG 
width is visibly poorer, and the log likelihood is worse than the optimum fit 
by 302, an enormous amount. An adequate fit requires a b\ width of at least 
180 MeV to reproduce Fig. 2(b). DM2 found a similar result [1]. However, 
changing the width to 142 MeV has almost no effect on the parameters fitted 
to the a amplitude. The reason for this is as follows. The o band crosses the 
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Fig. 2. The projection of the fit compared with data in Method I, using Eqns. 
(10)-(12) for the a parametrization. (a) and (b) show the projections of tttt and 
uj-k mass. Histograms show the fit, and the shaded region indicates the background 
estimated from sidebands. The dashed curve in (b) shows the fitted &i(1235) sig- 
nal (two charge combinations), (c) and (d) show the mass projections of the ++ 
and 2 ++ contributions to ir + Tr~ from the fit. In (c), the shaded area shows the a 
contribution alone, and the full histogram shows the coherent sum of the a and 
/o(980). 

two &i(1235) signals on the Dalitz plot. The fit effectively integrates over the 
whole hi band, and interferences with the a are affected little by the precise 
line-shape of the &i(1235). 

When using the parametrization of Eqn. (13) [8] for the a pole (Fit B), we 
find an optimum at M = 526 ± 15 MeV, r = 535 ± 50 MeV. This gives a 
pole position of M = (570 ± 7 (stat) ±19 (sys)) - i(27A ± 14 (stat) ±22 (sys)) 
MeV, in satisfactory agreement with Fit A. 

The Breit-Wigner amplitude of constant width for the a given in Eqn. (9) 
gives a very similar intensity distribution to that of Fig. 2. The mass and width 
optimize at M = 470 ±20 MeV, F = 613 ±60 MeV. Again the mass and width 
are strongly correlated; the errors cover these correlations. Table 1 summarizes 
the a pole positions obtained with Method I. In all cases, convergence of the 
fit is rapid, and the solution is unique. 



Angular distributions are shown in Fig. 3 for different slices of tnr mass. The 



BW Function 



Pole position (MeV) 



Eqns. (10)-(12) 542 ± 7 ± 15 ± 30 (extrap) - i(249 ± 15 ± 20 ± 30 (extrap)) 
Eqn. (13) 570 ± 7 ± 19 - i(274 ± 14 ± 22) 

Eqn. (9) 542 ± 7 ± 20 - i(269 ± 15 ± 25) 



Table 1 

Pole positions of the a for three fits from Method I. Here, the first error is statistical, 

the second is systematic and the last error (30 (extrap) MeV) accounts for the errors 

in the extrapolation to the pole, which can also be applied to other two fits in the 

table. 

third distribution, cosa^ departs significantly from isotropy. This effect was 
observed in the earlier DM2 data [1]. Up to M(tt7t) = 800 MeV, the departure 
from isotropy is due entirely to interference with the &i(1235). 



^5 500 

S 

63 



I I I I I I I I I I I I I I I I I 



m I r 



1000 



500 r 



.i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1, 



-I I— ■— i— i— I i i i I I I i I IT+J I - 



L r^^"^ J 500 N^«=^! 500 



100 



-10 1 

cosB„ 







I I I I I I I I I I I I I I I I I I I. 



Ill I I I IT 



I I I I I I I I I I I I I I I I I I I. 



-1 



Q In I 



-1 



cos a 



7U 



cos B, 

~7t 



a 

« 500 

s 

63 



,l l l l l l l l l I l l l l l l l. 



1000 - 



[a) 300 to 500 MeV 

1000 i_i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ij 1000 i ■ iiiiiiiiiiiiiiiiii_ 



I I I I I I I 1 1 I I I I I I I I I I 



T 



100 



J^X 



■10 1 

cosQ„ 




cos a. 



(b) 500 to 700 MeV 



Fig. 3. Angular distributions for angles %■> @u, Ott and f3 n in Method I when Eqns. 
(10)-(12) are used for the a parametrization. Here, x is t ne azimuthal angle be- 
tween the production plane of J/ip — > u>X and the decay plane X —> -kit, d u is the 
production angle of the u> in the J/tp rest frame, a f is the decay angle of the 7r + in 
the rest frame of X, taken with respect to the direction of the recoil u>, and /?„■ is 
the angle of the ir + with respect to the direction of X in the rest frame of the lo. 
Histograms show the fit. Slices of -kit mass are (a) 300-500 MeV and (b) 500-700 
MeV. 

The changes in log likelihood after removing or adding some components in 
the fit are studied with coupling constants of all other contributions being 
reoptimised. Our definition of log likelihood is such that a change of 0.5 cor- 
responds to a one standard deviation change in one fitted parameter. Fig. 4 
shows the fits when some of these components are removed. Fig. 4(a) is the 
fit without / (980). The log likelihood is improved by 224 when / (980) is 



10 



added. Therefore, /o(980) is required. It contributes about 1.1% of the inten- 
sity. Errors quoted above for the a pole cover the changes when the / (980) 
is removed. A fit without /2(1565) is shown in Fig. 4(b); the log likelihood 
is worse by 210. This fit fails to reproduce the dip in the mr mass spectrum 
at 1560 MeV. Alternative fits using / (1500) fail to reproduce this feature. 
Therefore, a small but definite contribution from /2(1565) is required with a 
fitted mass optimising at 1540 MeV with a width fixed at 126 MeV. 

A 7r + 7r~ contribution in the mass range 2000-2250 MeV is needed. Spin 2 
gives the best fit, and a good fit can be obtained using the / 2 (2240), listed by 
the PDG under ^2(2300) and reported by Crystal Barrel with a width of 240 
MeV [16]. However, it lies right at the top of the available mass range, so this 
identification is not unique. The fit with /2(2240) removed is shown in Fig. 
4(c). 

A definite contribution is required from p(1450). Fig. 4(d) shows the fit omit- 
ting this amplitude; it makes the log likelihood value worse by 923. Note that 
there are contributions to this mass projection from two ujtt contributions; this 
is why removing the p(1450) affects also the high mass region. 

We have tried adding a further broad contribution to the 7T7T S-wave with mass 
1500-1600 MeV and a large width of order 800 MeV; such a contribution 
has been found by Anisovich et al. in fitting other data [18]. This gives no 
significant improvement. If /o(1500), /o(1710) or / (2100) are added, there is 
no optimum if their masses are scanned. Fitted contributions are 0.43% for 
w/ (1710), 1.1% for w/ (1500), and 0.36% for / (2100). We regard these as 
upper limits and omit them from the final fit. 



4.2 Method II 



The second independent analysis, Method II, described in Section 2, is now 
reported. Three parametrizations of the a pole - Eqn. (9), Eqn. (13), and Eqn. 
(14) are used. The background is fitted by a non- interfering 57r phase space. 
The width of the a particle is the width at its mass, i.e., r cr (m -). Masses and 
widths of the a particle are obtained from the optimisation. 

Fig. 5 shows the projection of the fit compared with data when using a Breit- 
Wigner amplitude of constant width (Eqn. (9)). The mass and width of the a 
pole optimise at M a = 446+^2° MeV and T a = 578±1±86 4 MeV - 

The masses, widths, and pole positions of the a are shown in Table 2 when 
different a parametrizations are used in Method II. The first errors are sta- 
tistical errors, and the second are systematic, which are determined from the 
variation for different treatments of the background (sideband subtraction or 
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Fig. 4. Method I: (a) Fit without /o(980), (b) Fit without / 2 (1565), (c) Fit without 
/ 2 (2240), and (d) Fit without p(1450). 

direct fit), the changes in the solution when adding or removing small compo- 
nents, as well as the differences when changing the background level. The final 
global fit to the angular distributions of M n + n - < 1.5 GeV/c 2 mass region is 
shown in Fig. 6. The fit agrees well with the data. 



BW Function 


Mass (MeV) 


width (MeV) 


Pole Postion (MeV) 


Eqn. (9) 


446±^ 3 f 


C7S+36+114 
J '°-23-86 


(512±Jg±« ) - i (252 +££> ) 


Eqn. (13) [8] 


530^+f 


448±f 7 ^ 9 


(558±#±£) - i (231 tUttl) 


Eqn. (14) [17] 


7CO+10+76 
' oz -10-77 


984±|S±lt| 


(521±H±tf ) - i (237 i?i|) 



Table 2 

Masses, widths and pole positions of the a particle for Method II. The first errors 

are statistical, and the second are systematic. 

Both the a and the /o(980) are obviously needed in the data. Omitting the a 
makes the log likelihood worse by 5238. If the /o(980) is removed, the log likeli- 
hood is worse by 202. The largest resonance in this channel is the ^2(1270) with 



a mass and width of M 



/ 2 (1270) 



1268 ± 4 MeV and T 



/ 2 (1270) 



180 ± 10 MeV, 

respectively. Another important resonance in this channel is the &i(1235). 
The vertical and horizontal bands in the Dalitz plot (Fig. Id) correspond to 
J/ip — > 6i(1235) ± tt =f . In the 7r + 7r _ mass spectrum, the 6i(1235)tt contribu- 
tion extends below the /2(1270). They interfer with the a-particle. The mass 
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Fig. 5. The projection of the fit in Method II when Eqn. (9) is used as the a 
parametrization compared with data,, (a.) The shaded histogram is all ++ con- 
tributions including the a, the /o(980), and their inteference. (b.) The shaded his- 
togram is all 2 ++ contributions, (c.) The shaded histogram shows the contribution 
from the 6i(1235). (d.) The lower shaded area represents the background, where the 
curve shows the fitted phase space background and the histogram is the background 
estimated from u> sidebands. 

and width of the &i(1235) are M = 1231 ± 12 MeV and T = 244 ± 24 MeV, 
respectively. As in Method I, the width of the &i(1235) is much larger than 
that of the PDG. However, fixing the width of the &i(1235) to the PDG value 
does not change the parameters of the a appreciably. The contribution of the 
background is shown in Fig. 5(d). 



5 Summary 



In summary, the essential features to emerge from these data are the existence 
of the a and the determination of the a pole position. Two independent anal- 
yses are performed, and different parametrizations of the a pole are applied. 
The mass and width of the a are different when using different a parametriza- 
tions. However, the pole position of the a is stable. Different analysis methods 
and different parametrizations of the a amplitude give consistent results for 
the a pole. From a simple mean of the six analyses, the pole position of the a 
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Fig. 6. The angular distributions in the mass region M w + W - < 1.5 GeV in Method 
II when Eqn. (9) is used for the a parametrization. Here, <f>\ is the azimuthal angle 
and 9\ is the polar angle of X in J/ijj rest frame for J/ijj — > uX, and 02 and 62 are 
the corresponding angles of the ir + in the X rest frame. The histogram is the fit, 
and the crosses are data. 



is determined to be (541 ± 39 - % (252 ± 42)) MeV. Here, the errors from the 
constant width Breit-Wigner parametrisation in Method II are chosen, which 
are larger than the errors of the fitted results in Method I. The systematic 
errors dominate. 
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